Localization and derealization in dirty superconducting wires 
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We present Fokker-Planck equations that describe transport of heat and spin in dirty unconven- 
tional superconducting quantum wires. Four symmetry classes are distinguished, depending on the 
presence or absence of time-reversal and spin rotation invariance. In the absence of spin-rotation 
symmetry, heat transport is anomalous in that the mean conductance decays like l/\/X instead 
of exponentially fast for large enough length L of the wire. The Fokker-Planck equations in the 
presence of time-reversal symmetry are solved exactly and the mean conductance for quasiparticle 
transport is calculated for the crossover from the diffusive to the localized regime. 

PACS numbers: 72.15.Rn, 73.20.Fz, 73.23.-b, 74.25.Fy. 



The discovery of the d-wave nature of the order pa- 
rameter in high T c materials has renewed interest in un- 
conventional superconductors with low energy quasipar- 
ticles near the Fermi energy Ep. An important question 
is how disorder affects the quasiparticle dynamics and 
the corresponding low-temperature properties of the su- 
perconductor. In Ref. Q it is predicted that, on energy 
scales \e — sp\ less than the inverse mean free time and 
on length scales beyond the mean free path, weak impu- 
rity scattering leads to a finite density of states (DoS) 
and to a diffusive dynamics of quasiparticles. In normal 
metals, it has been known for a long time that quantum 
interference imposes corrections to this picture, in the 
form of weak localization, and eventually, for dimensions 
< 2, exponential (Anderson) localization. The analogous 
question for low energy quasiparticles in unconventional 
superconductors has been considered only recently |^| . 

The crucial distinction between quasiparticles in a su- 
perconductor and in a normal metal, is that the for- 
mer are described by a Hamiltonian of Bogoliubov-de 
Gennes (BdG) type. Such a Hamiltonian has an addi- 
tional particle-hole grading, accompanied by a discrete 
particle-hole symmetry, which is absent in the Hamilto- 
nian for (electron-like) quasiparticles in a normal metal. 
Symmetry plays a crucial role in the problem of Ander- 
son localization. A classification of the symmetry classes 
for BdG Hamiltonians, depending on the presence or ab- 
sence of time-reversal (TR) and spin-rotation (SR) sym- 
metry, has been given by Altland and Zirnbauer || . The 
four possibilities are denoted C, CI, D, and Dili, see ta- 
ble 0. Ref. [§ addressed the "zero-dimensional" (0D) 
case of chaotic quantum dots with superconducting leads. 
The higher dimensional realizations of the BdG symme- 
try classes, relevant for the question of localization, were 
studied in Refs. mainly by field-theoretical methods 

involving construction and analysis of non-linear sigma 
models with appropriate symmetries. 

In this Letter we study localization in the BdG sym- 



metry classes for the geometry of a quantum wire, i.e., in 
quasi-one-dimension (quasi-lD). For this purpose, we use 
the Fokker-Planck (FP) approach |^,^0), which is comple- 
mentary to the non-linear sigma model of Refs. 
Using the classification scheme of Ref. j8| , we obtain four 
FP equations that control quasiparticle transport at the 
Fermi level in a dirty superconducting wire. Our findings 
are remarkable: While for classes C and CI the mean 
and typical values of the quasiparticle conductance g de- 
cay exponentially with the length L of the wire for large 
L, the situation in classes D and Dili is quite different. 
There the mean (g) decays only algebraically to zero for 
large L and lag is not self-averaging, indicating a very 
broad distribution of the conductance and the absence of 
the exponential localization of the quasiparticle states at 
ep. (The absence of exponential localization for class D 
has been announced independently in Ref. 0.) 

It should be stressed that the BdG Hamiltonians do 
not conserve charge. Instead, the conserved densities are 
those of the energy (in all four classes) and spin (when 
the SR symmetry is present). Thus, the transport prop- 
erties (the conductance g) studied in this Letter refer to 
transport of heat and spin. 

We now proceed with a detailed statement of our re- 
sults and their derivation. The model that we consider is 
that of a disordered quantum wire, with a Hamiltonian 
of the BdG form. We distinguish gradings corresponding 
to spin up/down, particle/hole, left/right movers. Denot- 
ing these with Pauli matrices a, 7, and r, respectively, 
we write our model Hamiltonian as 



H = JC + V, 



K = iv F d x a <x> 70 <E> t 3 ® 1 N , (1) 



where do is the 2x2 unit matrix in the spin grading etc. 
The kinetic energy K, describes the propagation of right 
and left moving quasiparticles in N channels at the Fermi 
level. The "potential" V(x) is an 8N x 8iV matrix that 
accounts both for the presence of disorder and of super- 
conducting correlations. In particle/hole (7) grading it 
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TABLE I. Definition of the symmetry classes for a dirty 
superconducting quantum wire, with respect to the presence 
of spin-rotation (SR) and time-reversal (TR) symmetry. The 
table lists the Lie group C of the transfer matrix M, the factor 
group Q of angular degrees of freedom of M, the multiplicities 
of the roots of the symmetric space C/Q, and the degeneracies 
of the radial coordinates Xj of the transfer matrix A4. 
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-A* 



A 



(2) 



where v (A) is a hermitian (antisymmetric) 47V x 47V 
matrix, representing the impurity potential (supercon- 
ducting order parameter). The form (|2|) of the poten- 
tial V ensures that the Hamiltonian Ti obeys particle- 
hole symmetry, Ti = —jiH T "fi M. In addition, Ti (and 
hence V) may obey TR invariance Ti = TTi*T~ l , with 
T ' = it\® (72, and/or SR invariance Ti = —j2Ti T j2- 

Spatial fluctuations of the order parameter A and the 
potential v are taken into account by assuming that V 
is a Gaussian random variable with vanishing mean, i.e. 
that its probability functional P[V] is of the form 



P[V] oc exp 



dxtrV 2 {x) 



(3) 



where £ is the mean free path, 7 is a numerical constant to 
be defined below, and c = 1 (2) for class C/D (CI/DIII). 
The transport properties of the Hamiltonian Ti describe 
the transport of spin and heat by quasiparticles in a dis- 
ordered superconducting quantum wire. 

Before we continue with the analysis of our model (Q- 
^), some remarks about its validity and relevance are 
in place. One key property of the model is that, apart 
from corrections at very low energies due to quasiparti- 
cle localization M or the appearance of a critical state, 
the DoS of the Hamiltonian Ti near the Fermi level £f 
is nonzero and finite. This is related to the fact that 
the statistical average of the order parameter A is zero 
in our model, cf. Eq. ([|). For a dirty superconductor, 
such behavior is plausible if the order parameter is un- 
conventional, as in d-wave superconductors, or when it 
breaks TR symmetry, as is believed to be the case for, 



e.g., the ruthenates |11|, vortex lines in a (conventional) 
superconductor |3| , or a normal metal wire with magnetic 
impurities that is weakly connected to a superconduct- 
ing substrate. In all these cases the disorder leads to the 
existence of low-energy quasiparticle states plp[] . The 



Hamiltonian ([j]) then describes diffusion and localization 
of these "disorder-facilitated" quasiparticles. An alto- 
gether different scenario is that of a wire made out of an 
unconventional superconductor with very weak disorder. 
If boundary conditions are suitably chosen, one or sev- 
eral propagating modes can exist at £p, whose localiza- 
tion properties are described by Eq. ([!]). In any case, one 
should view Ti as an effective or coarse grained Hamilto- 
nian, whose validity is restricted to length scales beyond 
the microscopic mean free path £. It is universal in the 
sense that its form is determined solely by the symmetry, 
and the distribution (J|) provides for the existence of the 
diffusive regime with a finite DoS at the proper energy 
scale. (Note that the restrictions to the validity of our 
model are not different from those of related field theo- 
retic descriptions appearing in the literature H^HR.) 

We describe transport properties of the model (jl]-|§) 
through its 87V x 87V transfer matrix M. that encodes the 
x-dependence of an 87V-component quasiparticle wave- 
function tp satisfying the Schrodinger equation Tiip = eip 
at e = 0, ip(x + L) = M(x + L,x)t()(x). Formally, M is 
related to the Hamiltonian ([!]) as 



Ai(x + L; x) — T y exp 



dyr 3 V(y) 



(4) 



where T y denotes the path ordering operator for the 
y-integration along the wire. From Eq. (^) one finds 
that flux conservation (i.e., Hcrmiticity of Ti) and 
particle-hole symmetry imply that M^t^M. = r 3 and 

71. M71 = M.* , respectively. Further, TR invariance re- 
quires TAAT^ 1 = Ai*, while SR invariance is obeyed if 

72. M72 = A4*. The transfer matrix A4 obeys the multi- 
plicative rule M(z,x) — M(z 7 y)M(y,x) for x < y < z 
and hence is an element of a certain Lie group C. The 
appropriate Lie groups for the four symmetry classes are 
listed in Table ||. We note that the actual transfer ma- 
trix group is an 87V-dimensional representation of the Lie 
group £, where C also allows a lower dimensional (irre- 
ducible) representation for the classes C, CI, and Dili 
|l2| . Elements of C are conveniently parameterized in 
terms of their polar decomposition, which, in an irre- 
ducible representation, takes the form 
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Here V, G 0(47V) [Sp(7V)] for classes D/DIII [C/CI], for 
all i = 1,2,3,4, and X is a diagonal matrix with posi- 
tive entries Xj. (By Kramers' degeneracy, the elements 
of X occur in pairs in class C.) The Xj serve as radial 
coordinates on the Lie Group C. One verifies that the 
eigenvalues of the true 87V x 87V transfer matrix M.Jvft 
occur in d-fold degenerate inverse pairs exp(±2xj), where 
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the degeneracy d is listed in Table Q. Hence the number 
of independent Xj 's is AN/ d. Finally, we note that the Xj 
are related to the conductance g through |1C| 



4N/d 

d coslT 
3=1 



(5) 



Our aim is to find the probability distribution of the 
Xj for a transfer matrix corresponding to the model (|l|- 
||). Increasing the length L of the wire by a small 
increment SL amounts to multiplication of its trans- 
fer matrix M(L) = M(x + L,x) by a transfer matrix 
M' = M(x + L + SL,x + L). Since M' is close to the 
unit matrix, random, and statistically independent from 
M(L), we find that as a function of L, M(L) performs 
a random trajectory on its Lie group C. Actually, we do 
not need to know the full trajectory on C if we are only 
interested in the conductance g. It is sufficient to know 
the trajectory of the radial coordinates Xj of M(L) after 
dividing out a maximal compact subgroup Q of C corre- 
sponding to the angular degrees freedom of M that leave 
the product MM' invariant, or, in other words, to know 
the trajectory of the Xj in the symmetric space C/Q [ PI . 
The subgroups Q are listed in Table |. Starting from the 
microscopic model (HH), one can show that the trajec- 
tory obeyed by the Xj IS cL Brownian motion on the coset 
space C/Q described by the joint probability distribution 
P(xx, . . . , x iN f d ; L). The L-evolution of P is described 
by a FP equation, which follows either from a direct cal- 
culation starting from Eq. (Q), or from the general theory 
of symmetric spaces pa] . In both cases we find 



dL 2-f£ 



4N/d 

- y 



3=1 



dxj 



1 I Aj-xp 

dxj 



(6) 



4N/d 4N/d 

J=Y[ |sinh2a; ; ,-| mi ]J Y[ I sinh( Xj ± x k )\ m °, 

3=1 k>j ± 

where the numbers mi and m are the long and ordinary 
root multiplicities in the symmetric spaces C/Q, see Table 
| and 7 = (4Nm /d) + l-m + mi. The FP equation (|) 
is supplemented with the boundary condition dP/dxj = 
(P/ J)dJ/dxj at Xj = 0. The initial condition M = 1 for 
L = corresponds to P{x\, . . . ,x±N/d, 0) = FJj 5(xj). 

The FP equation (^J) is the fundamental equation that 
governs quasiparticle transport and localization in quan- 
tum wires of the symmetry classes C, CI, D, and DHL 

In the localized regime L » N£, typically all Xj 
and their spacings are much bigger than unity, and the 
conductance is governed by the smallest coordinate x±. 
That coordinate has a Gaussian distribution, with mean 
miL/j£ and variance L/j£. For classes C and CI this 
implies that g is exponentially small, with 



(In g) 



2m;L 

7 £ ' 



1 4L 
var In 5 = — , 



(7) 



with mi — 2 for class CI and mi = 3 for class C. Expo- 
nential localization for class C in quasi-lD was previously 
obtained by Bundschuh et al. 0], using the non-linear 
sigma model. For class D and Dili however, to/ = 0, so 
that there is no exponential localization. Instead, g has a 
very broad distribution (broader than log- normal), with 
an algebraic decay of the mean and the variance and an 
L 1 / 2 -dependence of Ing, 




2d 



var g 



(3), 



4(tt - 2)L 
var lng= — -—. (8) 
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Hence in classes D and Dili, quasiparticle states are not 
localized at the Fermi level. Since they are neither truly 
extended (typically g <c I in class D and Dili), we la- 
bel them critical, following terminology from the case of 
quantum wires with off-diagonal disorder, where similar 
behavior is found at the center of the band Jl4],[l5| . 

The effect of disorder is much less pronounced in the 
diffusive regime I <C L -C N£. Here, the conductance 
has only small fluctuations around its mean. Following 
the method of moments jl(J, we find (g) from the FP 
equation by construction of evolution equations for the 
moments of g a = d^2^ cosh _2a X j, a — 1, 2, . . ., 



j£ d (g a ) _ mo 
d 



a-l 



a dL 



— — 2_^\9a- n g n ) J- 2 ,\9a- n+lgn/ 

n— 1 n— 1 

+ (am a - 2a - I + mi)(g a+1 ) 

+ (2a - am a + m a - 2mi)(g a ). (9) 

In the diffusive regime one may replace the average of a 
product by the product of the averages, and hence one 
finds for £ < L < N£ 



(.9) = 



AN£ d(m Q - 2mi) 



L + £ 



3m D 



0{l/L,L/N£). (10) 



The first term in Eq. (|io| ) is the Drude conductance, while 
the second term is the first quantum interference correc- 
tion to the average conductance. For classes C, CI, D, 
and Dili it takes the values -2/3, -4/3, 1/3, and 2/3, 
respectively. The weak localization correction for class 
C was obtained earlier in Ref. For the classes D and 
Dili, the correction is positive, i.e., quantum interference 
enhances the conductance relative to the classical Drude- 
like leading behavior (see also Ref. ||). This is similar 
to the phenomenon of anti-localization in the standard 
symplectic symmetry class, though, as pointed out by 
Bocquet et al. 0, here it is a precursor of the break- 
down of exponential localization, while in the standard 
symplectic class localization takes over in higher order 
quantum corrections. 

In the presence of TR symmetry, Eq. (||) is soluble. 
This is in contrast to the case of the FP equations for 



3 



the standard and chiral symmetry classes, where only 
the case of broken TR symmetry was exactly solvable 
The first step is a map of Eq. (|) onto 
a Schrodinger equation in imaginary time for the wave 
function ^({xj}; s) — exp [— | In J({xj})] P({xj};s) for 
47V/ d fermions in one dimension with coordinates on the 
half- line x > 0. They interact through a two-body poten- 
tial proportional to m — 2 in the presence of a one-body 
potential proportional to (to; — 2)m;. Hence, for classes 
Dili and CI these fermions are free and they only differ 
by the boundary condition obeyed by their wave func- 
tions * at the origin. We thus find the solutions 



/2L 



P oc Y[( Xj sinh2a; 3 ) m! / 2 e- 7 ^ 

3 

x J^J {x\ — xj) (sinh 2 xu — sinh 2 Xj) . (11) 



Using the method of bi-orthogonal functions |ll| it is 
then possible to calculate the average conductance (g) 
for all TV and L. Here we report the result for the limit 
of large TV, leaving the results for finite TV for a future 
publication, 



(</} 



s 3 



n— 1 

where s = Lj (4N£). Note the agreement with Eq. ( |l0| ) in 
the diffusive regime s <C 1. In the localized regime s»l, 
Eq. ( |l2] ) may be resummed, and the asymptotic result (|^) 
is reproduced for class Dili, while for class CI one finds 
(g) = 8(7ts) -1 / 2 exp(— 4s). Quite remarkably, the exact 
results (|l2l ) for (g) in classes CI and Dili are related 
to the average conductance (g) c h m the chiral unitary 
ensemble for odd channel number M] as (g(s/2))ci + 
2(<7(s))diii = 4( 5 (s)) ch . 

The absence of localization in wires of classes D and 
Dili may have important implications for higher dimen- 
sions, provided our results can be extended beyond ID, 
and provided they are not restricted to the regime of weak 
disorder. With respect to the latter restriction, we can 
point to the close formal similarity of the derealization 
for the FP equations of class D /Dili and the correspond- 
ing FP equation for the chiral symmetry classes with odd 
TV, where it is understood that the absence of localization 
holds both for weak and strong disorder |l4|,[l7j . Thus ar- 
guing that quasiparticle states at the Fermi level remain 
delocalized for arbitrary disorder strength and dimen- 
sionality in the D-classes, our result suggests a possible 
resolution of a controversy in the literature surrounding 
2D disordered superconductors of class D j6jQ] . While all 
Refs. ||f7j assumed existence of two localized phases, dis- 
tinguished by the quantized value of the Hall conductiv- 
ity cr xy , and a metallic phase, the proposed global phase 



n 2 /2s_ 



CI, 



Dili, (12) 



diagrams and transitions between the phases differ con- 
siderably. We suggest that the solution might simply lie 
in the absence of localized phases for classes D and Dili 
in any dimension > 1. 

In conclusion, we considered quasiparticle transport 
and localization in disordered quasi- ID superconduct- 
ing wires at the Fermi level for the four Bogoliubov-de 
Gennes symmetry classes C, CI, D, and DHL We ob- 
tained and solved the Fokker-Planck equations for the 
probability of the radial coordinates of the transfer ma- 
trix. While quasiparticle states are localized in classes 
C/CI, localization is absent if spin-rotation symmetry is 
broken (classes D/DIII). 
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